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The branching fraction for the decay t^Vt, with —> e'^Vei'T, is measured using a data 

sample corresponding to an integrated luminosity of 427 fb~^ collected at center of mass energies 
neax 10.58 GeV with the BABAR detector at the PEP-II asymmetric-energy e"'"e~ collider at SLAG. 
In the process e"'"e~ — » cc — » D*s^DtagKX, the D*"*" meson is reconstructed as a missing particle, 
and the subsequent decay D*'^ — > D^7 yields an inclusive Dt data sample. Here Dtag refers to a 
fully reconstructed hadronic D decay, A" is a or K''\ and X stands for any number of charged 
or neutral pions. The decay Df KgK'^ is isolated also, and from ratio of event yields and known 
branching fractions, B{D'^ —^ t^Vt) = (4.5±0.5±0.4±0.3)% is determined. The pseudoscalar decay 
constant is extracted to be fo^ = (233±13±10±7) MeV, where the first uncertainty is statistical, 
the second is systematic, and the third results from the uncertainties on the external measurements 
used as input to the calculation. 



The Df meson can decay purely leptonically via an- clean probe of the pseudoscalar meson decay constant 
nihilation of the c and s quarks to a virtual boson fu^ , which describes the amplitude for the c and s quarks 
which decays to a lepton pair. These decays provide a to have zero spatial separation within the meson, a nec- 



essary condition for the annihilation to take place. In the 
Standard Model (SM), ignoring radiative processes, the 
total width is 



2jl 



M 



M2 
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where Mj~,+ and rui are the Df and lepton masses, re- 
spectively, Gp is the Fermi coupling constant, \Vcs \ is the 
magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element that characterizes the coupling of the 
weak charged current to the c and s quarks [1] . 

The Icptonic decay of the Df meson is hclicity- 
suppressed because it has zero spin, so that the final 
state neutrino and lepton must combine to form a spin- 

state. Consequently, the left-handed neutrino forces 
the anti-lepton to be left-handed, thus suppressing the 
decay rate by the factor m|/M^+. The net effect of 
hclicity and phase space factors results in large differ- 
ences in the leptonic branching fractions of the Df me- 
son. The branching fractions for Df decays to lu^, where 

1 = e"^, t"*", are roughly 2 x 10~^: 1 : 10 in propor- 
tion. The large branching fraction for the t"*" decay mode 
motivates the use of the decay sequence Df t^i't, 
r+ e^vJ^T in this analysis. The signal branching 
fraction B{Df — > r+i/r) relative to the well measured 
branching fraction B{D+ K^K+) = (1.49±0.09)% [2] 
is determined and used to extract the decay constant fjj^ . 

In the context of the SM, predictions for meson de- 
cay constants can be obtained from QCD lattice cal- 
culations [3 -8]. The most precise theoretical prediction 
for /tj^, which uses unquenched lattice QCD, is (241±3) 
MeV [5]. The most precise measurement of the branch- 
ing fraction for Df t'^Vt- (r"*" e'^Vf.V.r) yields 
B{D+ -> T+Vr)= (5.30±0.47±0.22)% [10] and the value 
fjj^ = (252.5±11.1±5.2) MeV. Decay constants of £> and 
B mesons enter into calculations of hadronic matrix ele- 
ments for several key processes and their theoretical pre- 
dictions. For instance the calculation of BB mixing re- 
quires knowledge of Jb- While leptonic decay of the B 
meson is heavily CKM suppressed, leptonic decay of D+ 
meson is CKM favored and thus resulting more precise 
measurements of fn^ can be used to validate the lat- 
tice QCD calculations that are applicable to S-meson 
decay. Several models involving physics beyond the SM 
can induce a difference between the theoretical predic- 
tion and the measured value. These include a two-Higgs 
doublet model [12] and a model incorporating two lepto- 
quarks [13]. It is important to have high precision deter- 
minations of 1 both from experiment and theory, in 
order to discover or constrain effects of physics beyond 
the SM. The Particle Data Group gives a world average 
of = (273±10) MeV [9] but this does not include the 
most recent results [10, 11]. 

Measuring the branching fraction for t^v-t re- 

quires knowledge of the total number of D'^ mesons in 



the parent analysis sample. Alternatively, the branching 
fraction can be measured relative to that for a Z?+ de- 
cay mode with well-known branching fraction, with the 
latter then used to obtain B{D^ t'^v^)-. this is the 
procedure followed in the present analysis. The decay 
mode £>+ — > ^7r+ has been used often as a normaliza- 
tion mode. However this is somewhat problematic, since 
determination of the branching fraction for this decay 
requires a Dalitz plot analysis of the D+ — > K^K~-k'^ 
process. A description of the Dalitz plot intensity distri- 
bution incorporates contributions from other quasi-two- 
body amplitudes such as K*{d,m)'^K+, K*{U'iQfK+ 
and /o(980)7r"'". The contributions of these other de- 
cay modes to the specific mass range used to define the 
(j) — > K^K~ rate have to be taken into account. These 
depend on the mass and width values of the resonances 
and their interference effects, as well as the mass reso- 
lution of the experiment. [2, 14]. For these reasons, the 
decay mode D+ i^g_ftr+ is chosen instead as refer- 
ence mode in the present analysis. Its branching frac- 
tion is quite well known, and the branching fraction for 
£>+ — > t^Vt can then be expressed as: 
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where Ns and e refer to the number of events and to- 
tal efficiency for the signal and the normalizing decay 
modes. The values of the branching fractions used for 
K° TT+TT- and r+ e+i^eVr are (69.20±0.05)% and 
(17.85±0.05)%, respectively [9]. 

This analysis uses an integrated luminosity of 427 
fb"""^ for e+e" collisions at center of mass (CM) ener- 
gies near 10.58 GeV, corresponding to the production 
of approximately 554 million cc events. The data were 
collected with the BaBMI detector at the SLAC PEP- 
II asymmetric-energy collider. The BaBAR detector is 
described in detail in Refs. [15, 16]. Charged-particle mo- 
menta are measured with a 5-layer, double-sided silicon 
vertex tracker (SVT) and a 40-layer drift chamber (DCH) 
embedded in the 1.5-T magnetic field. A calorimeter con- 
sisting of 6580 CsI(Tl) crystals is used to measure elec- 
tromagnetic energy. Charged pions and kaons are iden- 
tified by a ring imaging Cherenkov detector (DIRC) and 
by their specific ionization loss in the SVT and DCH. 
Muons penetrating the solenoid are detected in the in- 
strumented magnet flux return. 

The Df t'^Vt branching fraction measurement is 
carried out via the D*"^ production process e+e" ^ cc ^ 
D*+ DtagK'^'~ X , and the subsequent decay Z)*+ 
Df"/. Here, -Dtag is a fully reconstructed hadronic D 
meson decay, required to suppress the large background 
from non-charm continuum qq pair production; X repre- 
sents a set of any number of pions {n^ and tt^) produced 
in the cc fragmentation process, and JC^'" represents a 



single K'^ or K~ from cc fragmentation required to assure 
overall balance of strangeness in the event. The photon 
from the decay D*+ Df'j is referred to as the signal 
photon. 

Event selection begins with I^tag construction. Can- 
didates are reconstructed in the following modes: — > 
K+TT-{Tr'^), X+7r-7r-7r+(7r"), or i^-^Tr+Tr- (tt"), and 
^ K+TT-Tr-{Tr°), K°Tr-{Tr°), or i4:07r-7r-7r+. The 
probability for the geometric vertex fit of the TAG decay 
products must exceed 0.1%. The minimum required CM 
momentum of the X'tag candidate is 2.35 GeV/c. It is 
chosen near the kinematic limit for charm mesons aris- 
ing from B decays in order to eliminate the associated 
large combinatoric background. The mass of the £>tag 
candidate must lie in the range 1.7-2.1 GeV/c^. 

A single K~ or from cc fragmentation that does 
not have tracks in common with the Dtag combina- 
tion is found. Kaons are identified using information 
from the DCH and DIRC. A K'^ candidate is recon- 
structed through its decay to two charged pions which 
must originate from a common vertex. The dipion in- 
variant mass must be within 25 MeV/c^ of the nominal 
Kg mass value [9] , and the flight distance must be greater 
than three times its resolution. Neutral pions are recon- 
structed through their decay to two photons; the invari- 
ant mass of the photon pair must be within 10 MeV/c^ 
of the nominal 7r° mass value [9] . Any charged or neutral 
pion not associated with the Dtag or the fragmentation 
kaon is assigned to the fragmentation X candidate. 

A D*"*" candidate is reconstructed as the missing parti- 
cle with its four-momentum defined as, Pjj»+ = Pe+e- ~ 
{Pd^j^^+Pk%/k- +-Pa:), where the four-momenta (P) are 
from the initial state, I?taG) the fragmentation kaon and 
the fragmentation X, respectively. The mass of the 
candidate must be within 200 MeV/c^ of the nominal 
£>*+ mass value [9]. The production vertex of surviving 
candidates is then fitted using mass, energy and collision 
point constraints. 

In order to be consistent with the decay sequence 
— > T+i/x, r"*" — > e^VgVT^ a _D+ candidate is selected 
by requiring that there be a single e""" in the event. The 
6+ must have a minimum number of coordinates in the 
SVT and the DCH to ensure a good quality track. Sim- 
ilarly, the decay K'^K'^ is identified by requiring 
a single i4r° and pair. In addition, these candidates 
must not have tracks in common with Dtag or the frag- 
mentation kaon. 

The four-momentum of the Df candidate, for both the 
signal and the normalization mode, is defined as the recoil 
in the two-body decay — *■ Dfj, Pj^+ = Pj^*+ — 
P-y, where Pj is the four- momentum of the signal photon 
candidate, which must have CM energy greater than ICQ 
MeV. The resulting D+ candidate must have mass within 
200 McV/c^ of the nominal Df mass value [9]. 

Surviving Df candidates are further separated from 



background by requiring that the probability for the 
D*+ kinematic vertex fit exceed 0.1%, and that the 
CM momentum of the Df candidate exceed 3.0 GeV/c. 
The whole reconstruction procedure was evaluated us- 
ing GEANT4-based [18] Monte Carlo (MC) events gener- 
ated with EvtGen [19]. The generated MC samples for 
D+ T+Vr (r+ e+VeVr), D+ KlK+ and cc cor- 
respond to 14, 26 and 2 times the acquired data samples, 
respectively. 

After the final selection, the only background decay 
modes which contribute to the peak at the recoil 
mass (with the expected yields and shape determined 
from MC events weighted to 427 fb""'^) are D+ rjc^v^ 
(226 events), D+ ri'e+Ve (24 events), D+ (j)e+Ve 
(75 events) and D+ K^e+Ve (59 events). 

The yields for the signal and normalization mode arc 
determined from unbinned maximum-likelihood fits to 
the respective recoil D+ mass and extra energy (E'extra) 
distributions. As described earlier, the recoil D'^ four- 
momentum is defined as Pq+ = Pjj-'+ — P-y] E^xua. is re- 
constructed as the sum of the CM energy of all photons 
in the event with laboratory energy greater than 30 MeV 
which are not associated with any of the reconstructed 
charged-particle tracks or reconstructed neutral pions of 
the event. The signal photon is also excluded. The value 
of i^extra has been found to discriminate most effectively 
between signal and background events when it is required 
to be in the range 0-0.5 GeV. The distributions of -Ecxira 
for signal MC, background MC (events passing selection 
that do not include the signal) and data are shown in 
Fig. 1 . The difference between data and MC at i^oxtra = 
is due to MC under-estimation of beam-related back- 
grounds and of noise in the calorimeter. It has been veri- 
fied that the MC gives a a good description of the iSextra 
distribution in data for values of i^extra above 20 MeV. 
Correlations between the recoil Df mass and ^/^xtra were 
found to be negligible. Due to the discontinuity in the 
-Eextra distribution, the data are divided into two samples 
from which the results are determined using a simulta- 
neous unbinned maximum likelihood fit. For -Eoxtra = 0, 
only the recoil D'^ mass is used as a discriminating vari- 
able. For Eexua. > 0, EextTa and the recoil Df mass are 
used. The fit components are signal, peaking background 
and non-peaking background. The K^K^ mode 

is found to have no peaking background in MC, and thus 
no peaking background component is included in the fit. 

The signal recoil mass probability density function 
(PDF) consists of a bifurcated Gaussian function with 
a tail component (BFG) [20], plus a Novosibirsk fimc- 
tion [21]. The shape of the i^extra distribution for back- 
ground is taken from the data sidebands of the recoil Df 
mass distribution (Mpocoii < 1-95 GcV/c^ and MrccoH > 
2.0 GeV/c^). The remaining PDFs are empirical func- 
tions which describe the MC predictions and data. The 
background recoil mass PDF is a Novosibirsk function. 
The peaking background recoil mass PDF is a BFG plus a 




FIG 



1: Extra energy (-Ecxtra) in MC and data for _D+ — > r'^Vr 
(t'^ e^feVr)- The MC is normalized to the size of the data 
sample. The points represent data, the solid grey histogram 
is from background MC and the solid black histogram is from 
signal MC. The gap between and 20 MeV is due to the 
minimum energy requirement on photon candidates. 



Novosibirsk function. Second order polynomial functions 
are used for the signal and peaking background PDFs 
for E'cxtra- A Novosibirsk function is used for the back- 
ground -Ecxtra PDF. The parameters describing the shape 
of each PDF are obtained from fits to MC distributions. 
The systematic uncertainties introduced by this proce- 
dure are discussed below. Only the parameters speci- 
fying the number of signal and background events are 
allowed to vary in the fit. The number of peaking back- 
ground events is determined from the MC normalized to 
the data sample. Using MC pseudoexperiments, the fit- 
ting procedure is found to yield unbiased estimates of the 
signal yield. The fits to data are shown in Figs. 2-5, and 

yield Ng"^ =A4:8±36 events and Tvf °^^=333±28 events. 
Using Eq. (2) and the total efliciencies (e""" = 0.075%, 
^K°K+ ^ 0.044%) the branching fraction for t^Vt 
is measured to be (4.5±0.5)%, where the uncertainty is 
statistical only. 

The systematic uncertainties associated with the selec- 
tion criteria are evaluated by comparing the selection effi- 
ciencies for MC and data separately for each selection cri- 
terion. The selection efficiency is defined as Nn-i/Na\\ 
where N^-i is the number of events passing all selection 
criteria except the one being evaluated, and A'aii is the 
number of events passing all selection criteria. The un- 
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MC efficiencies for 



certainty is then defined as \1 — R 
where R is the ratio of data and 
each decay mode. The uncertainty associated with the 
probability for the kinematic vertex fit is 1.1%, and 
that with the CM momentum of the Z)+ meson is 2.7%. 
The uncertainties in the PDF distributions are evalu- 
ated by individually varying each PDF parameter by 



FIG. 2: Df candidate recoil mass for — > r^Vr (r"*" 
e^UeUT-) with -Boxtra = 0; the curve results from the fit de- 
scribed in the text. 
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FIG. 3: (a) Dt candidate recoil mass for Dt ~* t^Vt (t^ — » 

e^VeV-r) with Bextra > 0. (b) iJoxtra for i?cxtra > 0. The Solld 

curves result from the fit described in the text, and the dashed 
curves represent the signal contribution. 



one standard deviation and obtaining the correspond- 
ing variation in the fitted number of signal and back- 
ground events. To obtain the total uncertainty the indi- 
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assessed by varying the uncertainties for the individ- 
ual branching fractions, and refitting for the number 
of signal events. The peaking background modes and 
their branching fractions are: D+ rje^v^ (2.9±0.6)%, 
£1+ ^ r/'e+i/e (1.02±0.33)%, D+ cjie+Ve (2.36±0.26)% 
and D+ K^e+Ue (0.19±0.05)% [9]. Other sources of 
systematic uncertainty include tracking efficiency (0.34% 
per track) and e+ identification efficiency (0.82%). The 
uncertainties from tagging and fragmentation particles 
cancel in the ratio of the signal and reference modes. Ta- 
ble I summarizes the systematic uncertainty estimates on 
the branching fraction. 
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FIG. 4: D+ candidate recoil mass for D+ KIK+ {KI 
7r"'"7r~) with Bextra = 0; the curve results from the fit described 
in the text. 



TABLE I: Relative systematic uncertainty estimates on the 
branching fraction. 
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FIG. 5: (a) Dt candidate recoil mass for Dt K°sK+ 

{Kg TT + TT") with i^extra > 0. (b) EcxtT!^ for Bextra > 0. 

The solid curves result from the fit described in the text, and 
the dashed curves represent the signal contribution. 



vidual contributions are added in quadrature. The im- 
pact of the uncertainty in the number of peaking back- 
ground events for I?+ t'^Vt on the signal yield is 



Source 



Uncertainty (%) 



Event Selection 3.0 
Particle Identification 0.82 
Tracking 0.68 

Tv-r PDF Distribution +7.7 -4.7 

KIK PDF Distribution -1-4.9 -0.6 

Peaking Background -1-4.5 -4.3 



In conclusion, using an integrated luminosity of 427 
fb~^ collected with the BaBAR detector, the branch- 
ing fraction for the decay — > t'^Vt is measured to 
be (4.5±0.5±0.4±0.3)%, where the first uncertainty is 
statistical, the second systematic and the third from 
the uncertainties on the branching fractions for 
K",K+, K° TT+TT- and t+ ^ e+VeVr [9]. The de- 
cay constant is extracted using Eq. (1) and the val- 
ues of mT-(1776.84±0.17 MeV/c^), mu, (1968.49±0.34 
MeV/c^), TD^(0.500±0.007 ps) and by assuming \Vcs\ 
= \Vud\ = 0.97425±0.00022 [22]. The value obtained is 
fjy^ = (233±13±10±7) MeV, where the first uncertainty 
is statistical, the second is systematic, and the third is 
from the uncertainties on the external measured quanti- 
ties used in the calculation [9] . The results presented here 
agree to within one standard deviation with the most re- 
cent CLEO-c result for £)+ — > r+i^x [10, 11] and recent 
unquenched lattice QCD calculations of /d^ [5, 7, 8]. 
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